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1 Introduction 

The three dimensional Euler equations are evolution equations for the three 
velocity components u{x,t), 



coupled with a fourth equation, V ■ m = 0, expressing incompressibility. In 
this Eulerian formulation the velocity u and pressure p are recorded at fixed 
locations x G R'^. The velocities and pressure vanish at infinity or are pe- 
riodic. The pressure is determined using incompressibility. The equation is 
conservative and the total kinetic energy, / is a constant of motion. 

The Euler equations can be studied in terms of the vorticity {^)- The 
vorticity is a vector u = V x u corresponding to the anti-symmetric part of 
the gradient matrix Vw. It obeys a quadratic equation, whose nature is such 
that the magnitude of the vorticity may increase in time. If the amplification 
is not rapid enough then a well-known criterion (0) guarantees that no blow 
up can occur: if 



— + u -Vu + Vp = 



(1) 




sup \uj{x, t)\dt < oo 
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and the initial data are smooth, then the solution is smooth on the time 
interval [0,T]. The vorticity equation can be interpreted as the vanishing of 
a commutator 

[A,^^] = o (2) 

where 

d 
at 

is the material derivative and 

ft ^ CO -v. 

The characteristics of the first order differential operator Q are vortex lines; 
the characteristics of the material derivative Df are Lagrangian particle paths. 
The Lagrangian variables are the path maps a i— > X[a,t). The connection 
between the Lagrangian description and the Eulerian one is given by the 
relations 

u{x,t) = ——^-^^ ^ = 

In this paper we discuss a description of the Euler equations as a system of 
three coupled active vector equations. The description concerns Lagrangian 
quantities computed in Eulerian variables and thus combines the physical 
significance of the Lagrangian description with the analytical advantages of 
the Eulerian description. The description bears similarities to the Clebsch 
variable representation. The Clebsch variables are a pair of scalars, 6, ip that 
are constant on particle paths and can be used to re-construct the velocity 
via 

OXi OXi 

This interesting representation is somewhat restrictive: not all solutions can 
be represented in this manner. That is because the Clebsch variables impose 
special constraints on helicity. Helicity is the scalar product of velocity and 
vorticity h = u ■ u. Although h itself is not conserved on particle paths, the 
integrals 

/ h{x, t)dx — c 
Jt 

are constants of motion, for any vortex tube T. A vortex tube T is a time 
evolving region in space (not necessarily simply connected) whose boundary 
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is at each point parallel to the vorticity, u ■ u = where u is the normal to 
dT at X G dT. The constants c reflect the degree of topological complexity 
of the flow ( 1^ ) and in general are non-trivial but they vanish identically for 
flows that admit a Clebsch variables representation. Indeed, for such flows 
the helicity is the divergence of a field that is parallel to the vorticity h = 
— V - (na;). Topological properties of streamlines and vortex tubes are relevant 
to hydrodynamic stability ([|]) and turbulence (0, 0). The description of 
the flow that allows for arbitrary vortex structures is based on formula (p!5|) 
(0, H, 1], |T0|) that was used for numerical computations. Somewhat 



related Hamiltonian formulations have been introduced by several authors 

(ini, HI, 0, 111). 



2 Eulerian-Lagrangian Description 

The Lagrangian formulation of the Euler equations describes the flow in 
terms of a volume preserving diffeomorphism, the map a ^ X{a,t). The 
curve t X{a, t) is the Lagrangian path at label a and obeys Newton's law 

^^ = ^x(a,t). (3) 
The incompressibility condition for the map is 

detiVaX) = l. (4) 
The initial condition sets the labels at the initial time: 

X(a, 0) = a. 

The forces Fx in are 

Fx{a,t) = -(V.p)(X(a,t)) = - [(V„X(a, t))*]"' (Vap)(a, t) (5) 

with p{a, t) = p{X{a, t)) and where p is the Eulerian pressure. The notation 
M* means the transposed of the matrix M, (M*)~^ its inverse. Multiplying 
(D by (VaXy we obtain 

d^X(a,t) ,,,,,,,,, 

' {VaX{a,t)r = -(V.p)(a,t) (6) 



3 



or, on components 



d'^X^{a,t) dX^{a,t) _ dp{a, t) 
Pulling out a time derivative in the left-hand side we obtain 



(7) 



d_ 
di 



dX^{a,t)dX^{a,t) 



dt 



dai 



dq{a, t) 
dai 



where 

q{a,t) = p{a,t) 

We integrate (^) in time, fixing the label a: 

dX^{a,t) dX^{a,t) 



dX{a,t) 



dt 



dt 



where 



and 



n{a, t) 



M(o)(a) 



W(o)(a) - 

q{a, s)ds 

(9X(a,0) 
dt 



dn{a, t) 
dai 



(8) 
(9) 

(10) 

(11) 
(12) 



is the initial velocity. Note that n has dimensions of circulation or of kine- 
matic viscosity (length squared per time). The conservation of circulation 

follows directly from the form ([l0|). Let us consider 

A{x,t) = X-\x,t) 



(13) 



the "back-to-labels" map, and note that it forms a vector of active scalars 
(an active vector) 

dA 

D,A = — + u-VA = 0. (14) 
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Turning to (]I0|), multiplying by [{VaX{a,t))*] ^ and reading at a = A{x,t) 
we obtain the formula 

7/ ,\ / 1 / 4/ dn(x,t) 
u\x,t) = [ul,^{A{x,t))) - (15) 

where 

n{x,t) =h{A{x,t)) (16) 

The equation (|15D shows that the general Eulerian velocity can be written in 
a form that generalizes the Clebsch variable representation: 

u = {VAYB - Vn (17) 

where B = U(^q){A{x, t)) is also an active vector 

DtB = 0. (18) 

Conversely, and somewhat more generally, if one is given a pair of active 
vectors A = {A\x,t), ■ ■ ■ , {x,t)) and B = {B\x,t), ■ ■ ■ , B^^ {x,t)) of ar- 
bitrary dimension M, such that the active vector equations ([14|) and (|l^) 
hold and if u is given by 

M 

u{x, t) = Y^ B\x, t)V,A\x, t) - V,n (19) 

k=l 

with some function n, then it follows that u solves the Euler equations 

— + ti ■ Vm + Vtt = 
ot 

where ^ 

TT = DfU H — ImP 

Indeed, the only thing one needs is the kinematic commutation relation 

DtVJ = VM - {V.uyvj (20) 

that holds for any scalar function /. The kinematic commutation relation 
(|20D is a consequence of the chain rule, so it requires no assumption other 
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than smoothness. Differentiating ( |T9| ) and using the active vector equations 
|18D it follows that 



M 



n 



k=l 



-V,{Dtn) - (V,n)* 



M 



,fc=i 



3 The Active Vector Formulation 

The previous calculations can be summarized as follows: A function u{x,t) 
solves the incompressible Euler equations if and only if it can be represented 
in the form u = ua with 

u'Ax,t) = r'W.,t))?^^-?^ (21) 

and 

V ■ua = (22) 
where A{x, t) solves the active vector equation 

{dt + UA-V)A = 0, (23) 

with initial data 

A{x, 0) = X. 

The function (p represents the initial velocity and the function nA{x,t) is 
determined up to additive constants by the requirement of incompressibility, 
V ■ Ma = 0: 

A / N 9 \ ^^dA"'{x,t)\ 

The periodic boundary conditions are 

A{x + Lej.t) = A{x,t) + Lcj] nA(x + Lej,t) = nA{x,t) (24) 
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with Cj the standard basis in R^. In this case 

5A{x,t) = X - A{x,t), (25) 

nA{x,t), and UA{x,t) are periodic functions in each spatial direction. One 
may consider also the case of decay at infinity, requiring that 6a, ua and 
ua vanish sufficiently fast at infinity. The equation of state (|2l|, ^) can be 
written as 

UA = 'P {r {A{-,t)) VA™(-, t)} = P {{VAY 0(A)} (26) 

where 

P = 1 - VA" V ■ (27) 

is the Leray-Hodge projector (with appropriate boundary conditions) on di- 
vergence free functions. 

The Eulerian pressure is determined, up to additive constants by 

, , dnA(x,t) , , ^ , , 1, , ,,2 
V\x,t) = — h UA{x,t) ■ VnA{x,t) + -\uA{x,t)\ . 

The Jacobian obeys 

det (yA{x,t)) = 1. 

The vorticity 

loa{x, t) = V X Ua 
satisfies the Helmholtz equation 

DfuA = loa ■ Vma (28) 

and is given by the Cauchy formula 

ua{x, t) = [VA{x, t)]-^ aA{x, t)) (29) 

where C = V x is the initial vorticity. 

The advantage of an active vector formulation is that A has conserved 
distribution, that is, for any function $ 

j (^{A{x,t))dx = const; 

in particular ||y4(-, t) [[^^^^(rfj.) is constant in time. 



4 Local existence 



The proof of local existence of solutions to the Euler equations in the active 
vector formulation is relatively simple and the result can be stated econom- 
ically. The well-known local existence results in Lagrangian ([|15[) and Eu- 
lerian (0]) variables require more derivatives and thus use more restrictive 



function spaces. 

Theorem 1 Let cj) he a divergence free C^'^ periodic vector valued function of 
three variables. There exists a time interval [0, T] and a unique C([0, T]; C^'^) 
spatially periodic vector valued function 5{x^ t) such that 

A{x, t) = X + 6{x, t) 

solves the active vector formulation of the Euler equations, 

dA 

— + M- VA = 0, 

ot 

M = P{(VA(x,t))*(/)(A(x,t))} 

with initial datum A{x, 0) = x. 

The same result holds if one replaces periodic boundary conditions with 
decay at infinity. Differentiating the active vector equation (^) we obtain 
the equation obeyed by the gradients 

dxi dxj 
It is useful to denote 

P,i = 6,1 - d,A-'di (31) 

the matrix elements of the Leray-Hodge operator. Differentiating in the 
representation (p6|) and using the fundamental property 

P.^ = 

dxi 




we obtain 



CM); ^; ^ 

' dxi' dxi 



(32) 
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Recall that the function C is the curl of 0. This relation shows that the 
gradient of velocity can be expressed without use of second order derivatives 
of A and is the key to local existence: the equation ( ppj ) can be seen as a 
cubic quasi-local equation on characteristics. Let us make these ideas more 
precise. We will consider the periodic case first. We write C^'^, j = 0,1 
to denote the Holder spaces of real valued functions that are defined for all 
X G and are periodic with period L in each direction. We denote by 
ll/llo,;. the CO'^ norm: 

ll/llo,^ = sup + sup - f{y)\ f^^Vj (33) 

and by the C"^'^ norm: 

||/||i,,= ||/||o,, + i^||V/||o,, (34) 

where the notation | ■ ■ ■ | refers to modulus, Euclidean norm, and Euclidean 
norm for matrices, as appropriate. 

We break the solution of the problem in two parts, the map 6 —>■ u and 
the map u ^ §. We denote the first one W. 

W[5,(j)]{x,t) = P{{I + V5{x,t)y(t){x + S{x,t))} (35) 

This map is linear in but nonlinear in 6. 

Proposition 1 The map W[6,(j)] maps 

continuously. There exist constants C depending on /i alone so that 
||W[5,0]||o,M<C||0||o,,{l + ||V5||o,,}' 

and 

II VW^[5, 0]||o,;. < C|| V X 0||o,^ {1 + \\V6\\oJ' . 
hold for any 6 E {C^'''f, G {C^'^'f . 



9 



For the proof we note that W is made up from a number of operations. 
The first operation is the composition (l){x) ^ 4>{x + 5{x)). For a fixed 
5 G {C^'^Y map X I— > X + 5 is Lipschitz. Composition with a Lipschitz 
change of variables maps C^'^ into itself continuously (we say that it is a 
continuous endomorphism) . The joint continuity of [0, 5] (t){x + 5) in C^'^ 
follows naturally. The second operation is a sum of products of functions 
(a matrix applied to a vector). This is a continuous operation because the 
Holder spaces C-^'^, j = 0, 1 we chose are Banach algebras. The third and 
last operation is the linear operator P, which is bounded in Holder spaces. 
We need to consider also derivatives of W . We use the formula ( |32D and note 
that the expression for the gradient is made of similar operations as above 
and apply the same kind of reasoning. This finishes the proof. 

Time does not play any role in this proposition because the equation of 
state (5, (j)) W[5, 0] is time independent. The second half of the procedure 
does depend on time. Let us denote by 9 the map that associates to two 
continuous paths t i— >• and t i— *• 0(-,t) a new path t \—>- 9; the path 

t 9 = Q[6, (j)] is obtained by solving the partial differential equation 

an 

— + u-V9 + u = (36) 
ot 

where 

u = W[6i-,t),(Pi-,t)l 
periodic boundary conditions are imposed on 9 and zero initial data 

9{x,0) = 

are required. The Euler equation requires only the use of a time independent 
0, but allowing time dependent is very useful: one can thus treat more 
equations, in particular the Navier-Stokes equation. Let us consider the 
space 

of continuous (C^'^)'^ -valued paths defined on a time interval [0, T], endowed 
with the natural norm 

||^||i,P = sup||^(-,t)||i,^. 

t 
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We will consider also the weaker norm 

\\e\\o,r^sup\\e(.,t)\\o,,. 
t 

G is nonlinear in both arguments. 
Proposition 2 The map Q[5, (f)] maps 

Q : Vt y< Vt ^ Vt 

and is continuous when the topology of the source space Vt x Vt is the natural 
product C^'^ topology and the topology of the target space Vt is the weaker 
QO,^J. topology. Moreover, there exists a constant C depending on alone so 
that 

\\We{;t)\\o,^.< [J^ \\^u{;s)\\o,,ds^ {exp{C \\Vu{- , s)\\o,^ds}^ 
holds for each t <T with u = W[5, 0] and 6* = 9[(5, 0] . 

Proposition 2 states that the map © is bounded but not that it is continuous 
in the strong C^'^ topology. The proof follows naturally from the idea to use 
the classical method of characteristics and ODE Gronwall type arguments. 
Similar ideas are needed below in the the slightly more difficult proof of 
Proposition 3 and we will sketch them there and therefore we leave the details 
of the proof of Proposition 2 to the interested reader. 

In order to proceed let us take now a fixed (j), take a small number e > 
and associate to it the set 

I CVt 

defined by 

I - {S{x, t); 6{x, 0) = 0, \\V6{; t) ||o,^ < e, < T}. 

Combining the bounds in the two previous propositions one can choose, for 
fixed 0, a T small enough so that 

5^e[5,0] =S[5] 

maps 



11 



Inspecting the bounds it is clear that it is sufficient to require 

T||V X 0||o,^ < ce 

with an appropriate c depending on fi alone. Leaving 0, e and T fixed as 
above, the map S is Lipschitz in the weaker norm C^''^: 

Proposition 3 There exists a constant C , depending on ji alone, such that, 
for every ^2 ^ ^, the Lipschitz bound 

\\S[Si]-SM\o,v<C\\Si-52\\o,r 

holds. 

It is essential that Sj e T, so that they are smooth and their gradients are 
small, but nevertheless this is a nontrivial statement. An inequality of the 
type 

\\S[S,]-SM\o,v<C\\5,-52h,v 

is easier to obtain, but loses one derivative. This kind of loss of one deriva- 
tive is a well-known difficulty in general compressible hyperbolic conservation 
laws. The situation is complicated in addition by the fact that the consti- 
tutive law W depends on gradients. As we shall see, incompressibility saves 
one derivative. The heart of the matter is 

Proposition 4 Let (j) G (C^'^)"^ he fixed. There exists a constant depending 
on ji alone so that 

\\W[5^A] - WiS^AWW < C\\Si - S2\\o,M\\i,^ 
holds for any Sj e C^'^ with < 1. 

One could use the condition dj e C^'^ with < M but then C would 

depend on M also. 

Proof of Proposition 4. Denoting 

u=W[6,,(P]-W[62,<p], 

5 ^ 5i- 82, 

i^ix) = ^ {(j){x + 5i{x)) + (j){x + 52ix))) , 
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we write 

with 

and 

Now the bound 



v{x) = 4>{x + Si{x)) - (t){x + S2{x)), 
7 = ^{Si + S2) 



Mi = P{(V5)>} 
t^2-P{(I + V7)M 



||«2||o,^<C||5||o,^||(/)||l,;. 

is obtained in the same way as the bound in Proposition 1. (ActuaUy (p 
Lipschitz is enough here.) The dangerous term is Ui because it contains VS. 
But here we can "integrate by parts" and write 

ui = -P{(V^)*5} 

because of incompressibihty. The matrix is bounded in C"'^ and the 
bound follows again easily, as the bounds in Proposition 1. This ends the 
proof of Proposition 4. We draw the attention to the fact that the presence 
of the * (transpose) operation is essential for the "integration by parts" to 
be allowed. 

Returning to the proof of Proposition 3 we denote 6j = SSj, Uj — 
W{6j, (f)), u = Ui — U2, = 6i — 02 and write 

We consider the characteristics X{a,t) defined by 
dX U1 + U2, 

— = — ^ — (X,t), X((a,0) = a 

and note that in view of Proposition 1 and the assumption 5j G X, the char- 
acteristics are well defined for Q < t < T. their inverse A{x,t) = X-^{x,t) 
(the "back-to-labels" map) is defined too. Moreover, 

sup — <C 

t,a oa 
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and 

sup 

t,x 

holds with a constant C depending 

Solving by the method of characteristics we obtain 

e{x,t) = - F{X{A{x,t),s),s)ds. 
Jo 

Using Proposition 4 in conjunction with the bounds in Propositions 1 and 2 
we see that F{x, t) is bounded (uniformly in time) in C^'i^: 

sup\\Fi.,t)\\o,,<CU\U,M\o,v 

t 

Compositions with the uniformly Lipschitz X and A are harmless and we 
obtain the desired result 

ll^llo.P < C\\5\\o,v. 

This ends the proof of Proposition 3. The proof of Theorem 1 follows now 
using successive approximations. Starting with a first guess Si & I we define 
inductively 

5n+i = SSn e I. 

Proposition 3 implies that the sequence 5„ converges rapidly in the C'^''^ 
topology to a limit 6. Because X is convex it contains this weaker limit 
point, S E 1. Because S has the weak Lipschitz property of Proposition 3 
it follows that S5 — S. This actually means that A — x + 5{x, t) solves the 
active vector formulation of the Euler equations and that u — W[5, 0] solves 
the usual Eulcrian formulation. 

Now let us consider the case of decay at infinity. This case is instructive 
to look at this case because it illuminates the difference between (f), u, W on 
the one hand and x,5,Q on the other hand; the function spaces need to be 
modified in a natural fashion to accommodate this difference. The issue of 
decay at infinity is both a physical one - the total kinetic energy must be 
defined, and a mathematical one - P must be defined. But apart from this, 
the decay at infinity requirement does not hinder the proof in any respect. 



< C 



dA 
dx 

on u alone. Consider now the function 
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Theorem 2 Let ^ he a C^'^ velocity that is square integrahle 

j \(t){x)\^dx < oo 
and whose curl is integrahle to some power 1 < g < oo, 

y" |V X (l){x)\'^dx < oo. 

Then for e sufficiently small there exists a time interval [0, T] and a C^''* 
function S{x,t) such that 

sup||V5(-,t)||o,^<e 
t 

and such that x + S{x, t) solves the active vector formulation of the Euler 
equation. The velocity corresponding to this solution helongs to C^'^, is square 
integrahle and the vorticity is integrahle to power q. 

The proof foUows along the same hnes as above. Because 4> enters hnearly in 

the expression for W and because we control V5 uniformly, issues of decay at 
infinity of do not arise. In other words, the function space for velocities does 
not need to be a Banach algebra, rather a module over the Banach algebra 
of the S variables, which need not decay at infinity. 

5 The blow up issue 

Any solution of the Euler equation can be constructed using a sequence of 
near identity transformations. One starts out with 

<f) = UQ 

and solves for an interval of time <t <ti the active vector equation 

dA 

— +UA-VA = 
at 

A{x,0) = x. 
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At time t = ti one resets: 

(/) = Ui = UA{-,ti) 

and solves the system above again, for a new time interval ti < t < t2 and 
so one continues the solution. The local existence result guarantees that 

(tn+l - tn)||VM„||o,^ > C> 

and during this time the solution A{x, t) remains close to the identity in the 
sense that 6 = A — x obeys 

l|V5(-,t)||o,M<e 
with a prescribed e << 1. The formula ^ implies then that 

with a fixed K > 1. If the inequalities above would be sharp then, of course, 
the time steps would have to decrease exponentially and the procedure would 
diverge in finite time. It is possible that for certain initial data the bounds 
may be overly pessimistic and the solution may exist for a long time. But 
with the present knowledge, if one desires long-lived solutions for arbitrary 
three dimensional data then one needs to smooth either at the end of each 
step or during each time step. If one applies a smoothing procedure one 
evidently changes the problem and one introduces an artificial dissipation. 
There are many ways one could conceivable regularize the Euler equations. 
The physically correct energy dissipating equation is the Navier-Stokes equa- 
tion. Unfortunately it is not known in three dimensions if the Navier-Stokes 
equations have globally defined unique solutions that converge to solutions of 
the Euler equations. Even in two dimensions, where the existence of smooth 
solutions is known for both the Euler and Navier-Stokes equations, the sit- 
uation is not entirely trivial {^^, The two dimensional situation is 
characterized by the absence of vortex stretching. In the case of the three 
dimensional Euler equations the vorticity magnitude evolves according to the 
stretching equation 

Dti\uj\)=a\uj\. (37) 
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The stretching factor a is related to the vorticity magnitude through a prin- 
cipal value singular integral (|1T9[|): 



a{x,t) = P.V. I D{y,i{x,t),i{x + y,t))\uj{x + y,t)\y^^. (38) 



Here y is the unit vector in the direction of ^(x, t) = is the unit vector 
tangent to the vortex line passing through x at time t and D is a certain 
geometric factor. The geometric factor is a smooth function of three unit 
vectors, has zero average on the unit sphere, / DdSijj) = and vanishes 
pointwise when $,{x,t) = ±^{x + y,t). Because a has the same order of 
magnitude as dimensional reasoning suggests blow up of the type one 
encounters in the ordinary differential equation ^ = 



sup \uj{x, t)\ 



dt ' 
1 



T-t 



But if the vorticity direction ^ is smooth then a geometric depletion of a 
occurs; that means that a is of the order of magnitude of velocity times 
the magnitude of the spatial gradient of ^ (an inverse length scale, assumed 
to be finite). The two dimensional Euler equations correspond to the case 
^ = (0,0, 1) and a = identically. If 



/ sup \a{x, t)\dt < oo 



then no blow up can occur. This idea of geometric depletion of nonlinearity 
has been investigated theoretically and numerically for the Euler equations 
and for a quasi-geostrophic active scalar equation ([0, ||2^, p2|, 



[p^). In the Eulerian-Lagrangian formulation of the Euler equation the role 
played by smooth stratifications can be explained in the following manner. 
Consider functions w = of the form 

w^x, t) = (VA(x, t))X(y4(z, t)) (39) 

associated to arbitrary vectors i/j. Alternately, one might consider solutions 
of 

Dfw+{\/uAyw = (40) 



17 



with initial data -0. A particular example is provided by choosing ijj = (j), i.e. 
= wa 

WA = (VA)*0(A) (41) 

which obeys 

V X Wa = V X ua = uja- (42) 

Because the vorticity 

ua = V X ua (43) 

satisfies 

UA = (VA)-Y(A) (44) 

it follows that 



uja{x, t) ■ w^{x, t) = C(-4(x, t)) ■ ip{A{x, t)), 



(45) 



holds for any ip or, in other words 



Dt{ujA ■w) = 



(46) 

holds for any solution w{x,t) of (^Oj). Global regularity of a solution of the 
Euler equations would follow from (^Sj) if one could find a sufficient family 
of vectors ip. By a sufficient family for the initial velocity (j) and the time 
interval [0, T] we mean a family of vectors tfi such that there exists a non- 
negative function j{t) with f^ jdt < oo such that 

li^Ala;, t) I < 7(t) sup \uja{x, t) ■ w{x, t) \ 

holds for < t < T. A sufficient family for all two dimensional fiows is 
provided by just one ifj, ifj = (0,0,1) with 7 = 1. Generalizations would 
consist of situations in which one could find sufficient families that depend 
on the initial data and time and take locally the role played in 2D by the 
vertical direction. 

The blow up issue becomes, in terms of A, a question of formation of 
infinite gradients in conserved quantities. This is similar to the case of hy- 
perbolic conservation laws but with the significant difference that the under- 
lying characteristic fiow is volume-preserving: detCVA) = 1, the matrix VA 
is invert ible and 



((VA(x,t))- 



1 



rDet 



dA dA 



(47) 
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holds, where Cj = (Sjk) is the canonical basis in R^. Consider the Euler- 
Lagrange label differentiation 



"2 



dAkdAi\ d 



(48) 



From the commutation relation (^) and the A equation (|23| ) it follows that 

(49) 



holds for any j = 1,2,3. This commutation relation simply says that in 
Lagrangian coordinates, time and label derivatives commute. Note, from the 
formulas (^) and (g^) that 



e.puDet 



dA dA 
dxi ' dxi 



p 



(50) 



It is clear now that Df commutes with Qa = uja''^ because it is represented 



in terms of L^: 



= Q{A)Lf. 



Observe that, in view of the definition of the operators Lf, 

d 



it follows that 
on the other hand 



Lf = {{VAix,t))-' 
{VA{x,t))-' 



dxk 



Lf [xi]; 



Df{x, 



holds, so from the commutation relation (^91) we obtain 

Df{{VA{x,t))-') =Lf{v^A)- 



(51) 

(52) 

(53) 
(54) 

(55) 



This equation, which could have been derived also directly from (PPD, implies 
the vorticity equation because of (p9|): 

DfuA = Q{A)Lf{uA) = ^a{ua)- (56) 
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Because of the result in ([^) and (^), it is clear that the finiteness of 

^11 i\/Ai-,t)y^\\Loo(^ax)dt 
implies regularity. Or, using (^7[), we deduce that the finiteness of 

'''^l|vA(-,t)||ioo(,,)(it 

implies regularity. Let us introduce now the matrix 

C^ix, t; z) = {VA{x + t)),^ ((VA(x, t))-^) . . (57) 

and call it the Euler-Lagrange calibrator. The formula 

C^(x,t;^) = Lf"'*)(A(x + ^,t)) (58) 

shows that the calibrator measures the response of the Eulerian translation 
to an infinitesimal Lagrangian translation. Note that 

C^(x,t;0) =5,,. (59) 

The calibrator is a quotient of gradients at different locations and therefore 
locally spatially uniform, temporally arbitrary changes like dilations do not 
affect it. The vorticity equation ( |56[ ) can be expressed in terms of the Euler- 
Lagrange calibrator (||25[|): 

/ [ D(rr,^r n^\aJ7.\—\ — 

dxi 



AM = [ip-v-j fl(C.c-C.c:i)..,(i)|j| (60) 

Note that 



is bounded. It is therefore natural to conjecture that the smoothness of 
prevents finite time blow up for the Euler equations. This conjecture is true 
for the quasi-geostrophic active scalar. The interested reader is referred to 
(H) for details. 

The blow up question for the Euler equations remains open. Numerical 
calculations provide insight and hints, but the answer will have to be an- 
alytical. The considerations above point towards a possible incompressible 
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dispersive effect that hinders blow up: as the gradients of A become large 
the resulting rapid ([15D and non-uniform (|55D motion disperses the large gra- 
dients. This might cause instability of blow up or perhaps its suppression. 
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